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We examine the energy distribution in a three-dimensional model granular system consisting of 
N inelastic hard spheres contained in an open cylinder of radius R under the influence of gravity. 
Energy is supplied to the system by a vibrating base. We introduce spatially resolved, partial 
particle-particle "dissipations" Vp P (r) and T>p p (r) for directions parallel and perpendicular to the 
energy input, respectively. Energy balances show that the total (integrated) "dissipation" Dp P < 
while Dp p > 0. The energy supplied to the perpendicular directions is dissipated by particle-wall 
collisions. We further define (3 = — Dpp/D pp , which in the steady state represents the fraction of 
the power supplied by the vibrating base that is dissipated at the wall. We examine the dependence 
of P on the number of particles, the velocity of the vibrating base, the particle-particle restitution 
coefficient, and the particle- wall restitution coefficient. We also explore the influence of the system 
parameters on the spatially dependent partial dissipations. We observe that, unlike 2?| p (r), T>p P (r) 
changes sign. 
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INTRODUCTION 



In an isolated granular system of inelastic particles the 
energy-dissipating collisions rapidly lead to a total ces- 
sation of motion. If an energy source is present, such as 
a vibrating wall, a fiuidized steady state can be main- 
tained. The steady state is characterized by the balance 
of the energy input with the energy dissipated through 
particle-particle and particle-boundary collisions. Al- 
though there is a superficial resemblance, these non- 
equilibrium steady state granular systems have radically 
different properties from their equilibrium counterparts 
in systems of elastic particles. Examples include non- 
Maxwellian velocity distributions Q, ll H S H H 0, 
a lack of energy equipartition between rotational and 
translational degrees of freedom in a system of inelastic 
rough spheres Q, and a lack of equipartition of kinetic 
energy between the components of a granular mixture 
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Moreover, as has been observed in several studies, the 
breakdown of equipartition in granular systems extends 
to the translational degrees of freedom. In an experi- 
mental study of a three-dimensional vibrofluidized sys- 
tem Wildman et al. showed that the granular tem- 
perature in the direction parallel to the energy injection 
is larger than in the perpendicular directions. Morgado 
and Mucciolo [l^, using a variant of the Direct Simu- 
lation Monte Carlo method, demonstrated a similar ef- 
fect in a two-dimensional vibrated granular system. Pre- 
vost et al. [? ] showed that velocity correlations in vi- 
brated monolayers are different in the longitudinal and 
transverse directions. Barrat et al. [2(j studied a ver- 
tically shaken three-dimensional system and noted that 
the injected energy is transferred to the perpendicular di- 
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rections through particle-particle collisions. This results 
in apparent restitution coefficients that are greater than 
unity for some collisions, and lead Barrat et al. to pro- 
pose the Random Restitution Coefficient model in which 
the restitution coefficient of each collision is randomly 
sampled from a given distribution [2pl l2ll 122^ . 

Various theories of vibrated granular media have been 
proposed [II Hill III III |2Ej. Warr et al. [H| proposed 
expressions for the particle-particle dissipation and the 
power supplied to obtain a scaling relationship for the av- 
erage granular temperature. Their experimental results 
were not consistent with the theory, however. Kumaran 
|25j sought to improve the theory by using a Maxwellian 
velocity distribution. This had the effect of modifying 
the prefactors but not the exponents of the scaling rela- 
tions. McNamara and Luding |27| reviewed and extended 
the theories and compared their predictions to the results 
of event driven simulations. The comparison confirmed 
that the Kumaran scaling relationship applies only in di- 
lute systems. Sunthar and Kumaran [28j subsequently 
extended Kumaran's analysis to dense systems. More 
recently, Brey et al. [2^| used a hydrodynamic descrip- 
tion of an open vibrated system in two and three dimen- 
sions. They showed that the temperature profile exhibits 
a minimum as a function of height, and also obtained an 
accurate description of the energy dissipated. 

In this paper we explore the anisotropy of the en- 
ergy consumed or supplied by particle-particle collisions 
in a three-dimensional, vibrofluidized system of inelas- 
tic, hard spheres. More specifically, we resolve the to- 
tal energy supply or consumption into directions parallel 
and perpendicular to the energy source. Energy balances 
show that collisions consume energy in the parallel direc- 
tion and supply it to the perpendicular directions. We ex- 
amine the dependence of these quantities on the particle- 
particle and particle-wall restitution coefficients as well 
as the velocity of the vibrating base. We also examine 
the dependence of the local energy consumption/supply 
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FIG. 1: The model system of N particles of mass m and 
diameter d in an infinitely tall cylinder of radius R. The base 
is vibrated with an amplitude A and a frequency v. 



on position and the above system parameters. 



II. SIMULATION 



puter simulation conserve momentum, but dissipate en- 
ergy. Although several models have been proposed for 
inelastic collisions |2jJ, , we use the simplest approach, 
i.e. a constant restitution coefficient. This choice has 
produced results that are in good agreement with exper- 
iment 0, E3, E El E3- The restitution coefficient 
ranges between (a completely inelastic collision) and 
1 (an elastic collision). The post-collision velocities of 
colliding particles % and j are related to the pre-collision 
velocities (v; and Vj) by the following equations: 
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V; = V; - 



[(vi-vj) 



n n 



[(v ; - Vj ) -n)]n, 



(la) 



(lb) 



where c is the restitution coefficient and the unit vector n 
points from the center of particle i to the center of parti- 
cle j. Similarly, the equation for a particle- wall collision 
with a restitution coefficient of c,„ is 



v- = Vi - (1 + c w ) (vi ■ f)f, 



(2) 



where f is a unit vector that points from the center of 
particle i to the point of contact with the wall. Finally, 
the equation for a particle-base collision with a restitu- 
tion coefficient of q, is 



v- = Vi - (1 + c fc ) [(vi - v b ) • f] f , 

where Vb is the velocity of the base. 

Equations H an d 13 result in energy changes of 
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[( Vi - vj) -nY 



(3) 



(4) 



A. Simulation algorithms 

We use a three-dimensional, event-driven simulation 
of inelastic, hard spheres to model granular particles in 
a cylinder in the presence of gravity. Identical spheres 
with a diameter d and mass m are placed in an infinitely 
tall cylinder with a radius R (see Figure QJ. The simu- 
lation generates a sequence of particle-particle, particle- 
wall and particle-base collisions. The base vibrates with 
a symmetric, saw-tooth waveform with amplitude A and 
frequency v. This form was chosen for simplicity, and 
changing to another symmetric form, e.g. a sine wave, is 
not expected to have much influence on the steady state 
properties of the system [3fJ. The velocity of the base 
for this type of waveform is Vb = AAv. 

The simulations were initialized by randomly placing 
particles in the cylinder and assigning velocities selected 
from a Maxwell distribution. The particles were then al- 
lowed to evolve from the initial configuration to a steady 
state before data was collected for the average values pre- 
sented in this paper. 

The particle-particle collisions modeled in the com- 
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- (vi • r) 



(5) 



for particle-particle and particle- wall collisions, respec- 
tively. These energy changes are always negative for val- 
ues of c and c w less than 1. The energy change associated 
with a particle colliding with the base is 



m 
~2 



(l + c b ) 2 [( Vi -v b ) 



-TO (1 + C b ) (Vi • f ) [(Vi - V b ) • f] 



(6) 



which can be either positive or negative. Most collisions 
with the base impart energy to the particle, although if 
the collision occurs as the base is descending energy is 
lost. 



B. System Properties 

The simulated trajectories of the particles are used to 
calculate macroscopic properties of the system such as 
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TABLE I: Relations for converting between dimensionless 
(quantities with a star) and real (quantities without a star) 
quantities. 



and 



Conversion Equation 


Parameter 


7?* = R/d 


Radius 


H* = H/d 


Height 


A* = A/d 


Amplitude 


u*=u{d/g)^ 2 


Frequency 


v* = v/ (gd) 1/2 


Velocity 


T* = T / (mgd) 


Temperature 


E* = E/ (mgd) 


Energy 



the packing faction (77), the granular temperature (T a ) 
and the energy dissipation. 

To facilitate the comparison of the simulation results 
with different experimental systems we use dimension- 
less parameters. The characteristic length and mass are 
taken as the particle diameter d and mass m, respec- 
tively. A natural characteristic time is yj d/g where g 
is the gravitational constant. All numerical results pre- 
sented in this paper are dimensionless and can be con- 
verted to real units using the conversion factors presented 
in Table ID 

The packing fraction is defined as the ratio of the vol- 
ume occupied by particles to the sample volume: 
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where (v p ) is the average volume occupied by particles 
within a sample volume V, and (n) is the average number 
of particles in this volume. 

The granular temperature in the a direction is calcu- 
lated as 



T 



■■<«£> 



(8) 



where (v^) is the mean squared velocity in the a direc- 
tion. As shown in prev ious work, these are not equal and 
vary with position |l%ll&ll&l33| . The directional granu- 
lar temperatures for a system consisting of 2100 particles 
is shown in Figure 

To study the dissipation in the system we start with 
the overall energy balance: 

AE io +J2 AEf all +J2 AE^ + A^ + AK = (9) 



where n p is the number of particle-particle collisions, n w 
is the number of particle-wall collisions, and rib is the 
number of particle-base collisions. A$ and AK are the 
differences in the potential and kinetic energies, respec- 
tively, of the system at the start and the end of the sim- 
ulation: 
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(11) 



where e.g. z itini uai and v iiin iu a i are the initial values of 
the height and velocity of particle i, respectively. The 
final values are designated with the subscript final. 

A$ and AK are bounded and fluctuate around steady 
state values. The other sums in the total energy bal- 
ance, however, are unbounded as they increase linearly 
with the total number of collisions, n = n p + n w + rtt,. 
Therefore, for a sufficiently large number of collisions cor- 
responding to an elapsed time t Equation [5] simplifies to: 



D pp + D pw + P b = 0, 



where 



D 



i>ii 



(12) 



(13) 



is the dissipation due to particle-particle collisions, 

1 \ - 



D 



pw 



is the dissipation due to particle-wall collisions and 
1 



Ph 



AE" 



(14) 



(15) 



is the power supplied by the vibrating base. 

Motivated by the concept of directional granular tem- 
peratures for the system (e.g. Equation |HJ we separate 
the energy changes into the three degrees of freedom for 
translation. For example, the energy change in the z- 
direction due to a particle-particle collision is 



,|| _ 

ij 

rn 
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AE^ = — (1 + cY (v r • fi)' (A • 2) 



(1 + c) (v r • n) (n • z) (v r • z) , 



(16) 



where v r = Vi — vj is the relative velocity of the two 
particles and z is a unit vector pointing in the z-direction. 
A similar equation defines the energy change in the x- and 
y-directions, AE/j. The two energy changes are related 
by 



AEa = AEh + AEl < 



(17) 



We may then write two additional energy balances - one 
for the direction parallel to the energy input by the base 
(i.e. the z-direction), and another for the perpendicular 
directions (i.e. the x- and y-directions): 

£ AE l + £ + A$ + AK " = ( 18 ) 
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FIG. 2: The granular temperature (T a ) as a function of height and radial position in (A) the x and y directions and (B) the 
z direction. The contours designate a change in T a of 0.7. Simulation conditions: N = 2100, R = 14.5, A = 0.348, v = 1.13, 
c = 0.91, and c w = 0.68. 



and 

^A^ + ^A^+A^ = 0, (19) 

Tip n w 

where the superscripts || and _L denote the parallel and 
perpendicular directions, respectively. Since the particles 
are taken as smooth Ai?J = for collisions with the 
wall and AE^~ — for base collisions. Gravity acts in 
the parallel direction, hence the potential energy change, 
A$, appears only in Equation^] In the limit of a large 
number of collisions, Equations 1181 and 1191 become 

Dl + P b = (20) 

and 

Dp P + D pw = 0. (21) 

Since energy is never gained at the wall, D pw < 0, and 
on average input at the base, Pb > 0, it follows from 
Equations [2U| and that 

Dl < (22) 

and 

K > 0. (23) 



Like D pp , D pp is always negative (for inelastic collisions). 
Dpp, however, is positive, i.e. the collisions supply energy 
to the perpendicular directions. D^ p does not, therefore, 
represent dissipation of energy. The equality in Equation 
[53] applies only when the walls of the cylinder are elastic 
and AE^- = 0. Also since the total energy change due to 
a collision is either zero for an elastic collision or negative 
for c < 1 

D pp = Dl + L>4 < 0. (24) 

Thus, the amount of energy lost in the parallel direction 
is always greater than the energy gained in the perpen- 
dicular direction. 

Since we are interested in the amount of energy that 
is transferred to the perpendicular direction we find it 
convenient to define the fractional energy transfer 

J-s pP 

which by Equation [21 is bounded between zero and one. 
In the steady state it follows from Equations 1201 and 1211 

that [3 = —D pw /P}>. Thus, in the steady state (3 can 
be interpreted as the fraction of the power supplied by 
the vibrating base that is dissipated by particle-wall col- 
lisions. Similarly, the fraction of the power dissipated by 
particle-particle collisions is —D pp /Pb = 1 — (3. 
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We also investigate the spatial dependence of the 
particle-particle dissipation by introducing 

V vP^ z ) = \Y. AE ti ( 26 ) 

where the summation is over collisions that occur in the 
volume element at (r, z). A similar equation defines 

Dp P (r, z). The total system dissipation in the perpen- 
dicular direction is 

D^ = 2v / 25£,(r,z)rdrd* (27) 
Jo Jo 

A similar expression gives the relation between D pp and 
T>p P (r,z). In addition to T> pp (r, z) we also examine the 
reduced distributions 

V^ p (z) = 2n ( R V^{r lZ )rdr (28) 
Jo 

as well as 2?| p (z), which is defined similarly. 
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FIG. 3: Values of (3 for systems with 700 to 4200 particles 
and base velocities ranging from 0.786 to 15.7. Simulation 
conditions were R = 14.5, A = 0.348, c = 0.91, and c w = 0.68. 



III. RESULTS AND DISCUSSION 

We used the simulation to examine the dependence 
of (3 on the base velocity, the particle-particle restitution 
coefficient, and the particle- wall restitution coefficient for 
systems of varying particle number N . Additionally, we 
investigated the dependence of 2?| p (r, z) and D pp (r, z) on 
the packing fraction and temperature distribution in the 
system. Finally, we studied the trends in the radially 
averaged functions D pp (z) and T> pp (z) as the number of 
particles, the velocity of the base, the particle-particle 
restitution coefficient, and the particle-wall restitution 
coefficient were varied. Constant values of R — 14.5, A = 
0.348 and Cb = 0.88 were used in all simulations presented 
in this paper. 

A. Fractional energy transfer 

Figure [3] shows the effect of the base velocity on 
the transfer of energy between the parallel and perpen- 
dicular directions. Systems consisting of 700 to 4200 
particles were simulated with base velocities of vi, = 
0.786,1.57,3.14,7.86,15.7. f3 increases with increasing 
base velocity at constant N and decreases with increas- 
ing N at constant vt,. The former trend is due to the de- 
crease in density as the granular temperature increases, 
which leads to more particle-wall collisions. Thus, more 
of the energy input by the base is being dissipated at the 
wall. Increasing N at constant Uj increases the number of 
particle-particle collisions. The increase in the particle- 
particle collision rate reduces the fraction of the input 
power dissipated at the wall. 



It is interesting to note that j3 appears to approach 
an asymptotic limit as Vb increases, at least for 700 and 
1400 particle systems. We expect that the limit is de- 
pendent on the values of the particle-particle restitution 
coefficient c and the particle-wall restitution coefficient 
c w . The data in Figures 01 and [S] support this hypothesis. 

The fractional energy transfer was also calculated for 
systems of 700 to 4200 particles with particle-particle 
restitution coefficients ranging from 0.8 < c < 0.99. Fig- 
ure 01 shows the values of (3 as a function of c for each 
system size. As the particle-particle restitution coeffi- 
cient approaches unity the value of (3 increases. When 
c = 1 , according to Equation [21 (3 = 1 as all power 
must be dissipated at the wall. Figure 0] also shows that 
(3 decreases as the number of particles increases, as seen 
previously. 

We have devised phenomenological explanations for 
these trends. As the particle-particle restitution coeffi- 
cient increases, the granular temperature increases, and 
the particle bed expands leading to fewer particle-particle 
collisions and more particle-wall collisions as explained 
above. It appears that as the number of particles in- 
creases, (3 becomes independent of N. As the system be- 
comes denser, a greater fraction of the power is dissipated 
by particle-particle collisions. Eventually the power sup- 
plied is insufficient to sustain a fluidized system. The 
particles simply collapse and move as a whole on the 
vibrating base. This last situation cannot be studied di- 
rectly with the current version of the simulation. 

The effect of varying the particle- wall restitution coef- 
ficient at constant N and is shown in Figure[Sl As the 
restitution coefficient goes to one, (3 goes to zero. The 
elastic wall dissipates no energy, so the dissipation D p 
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FIG. 4: Values of 13 for systems with 700 to 4200 particles 
and particle-particle restitution coefficients ranging from 0.8 
to 0.99. Simulation conditions were R — 14.5, A = 0.348, 
v = 1.13, and c w = 0.68. 




FIG. 5: Values of /3 for systems with 700 to 4200 particles and 
particle-wall restitution coefficients ranging from 0.5 to 1.00. 
Simulation conditions were R = 14.5, A = 0.348, v — 1.13, 
and c = 0.91. 



goes to zero according to Equation [21 As c w decreases 
the wall is able to dissipate more energy. As in the other 
graphs, an increase in the number of particles is seen to 
decrease the value of (3. 



B. Local dissipation 

We investigated the dissipation as a function of posi- 
tion in systems containing 2100 particles (roughly three 
layers of particles resting on the base). We determined 



the packing fraction and the dissipation in the parallel 
and perpendicular directions as functions of height and 
radial position. 

FigureEK shows a contour plot of V p ■ (r, z) and Figure 

03 shows a contour plot of T>\ p (r, z). While the latter is 
everywhere negative, energy is both dissipated and gen- 
erated in the perpendicular direction, depending on the 
position. The dissipation shows radial and height de- 
pendencies for both the perpendicular and parallel direc- 
tions similar to those observed for the packing fraction 
shown in Figure [Bp. The extrema in the dissipation pro- 
files, however, occur closer to the base than the maximum 
in the packing fraction. The maximum in T> pp {r,z) oc- 
curs at the same location as the minimum in T> pp (r,z). 
It should be noted that the minimum in the parallel 
direction is larger in magnitude than the maximum in 
the perpendicular direction as required by Equation 1171 
The variations in T> pp (r, z) and T>\ p {r 1 z) as a function of 
height are different, as shown in Figure Specifically, 
D pp (r, z) is positive at small heights with a negative min- 
imum in the middle while V pp (r, z) is negative over the 
entire region, i.e., energy is never gained in the parallel 
direction from particle-particle collisions. It is interest- 
ing to note that the trends observed for V pp (r, z) and 

T> pp (r, z) appear to be the inverse of those seen in the 
temperature profiles presented in Figure [21 



C. Influence of the number of particles on local 
dissipation 

We varied the number of particles in the system from 
700 to 4200, corresponding to approximately 1 to 6 layers 
of particles at rest. The base velocity was kept constant 
at vj, = 1.57. 

The values of V^ p {z) and T>\ p (z) for each system are 
shown in Figure Q The dissipation in the perpendicu- 
lar direction exhibits a positive maximum for all the sys- 
tems studied. In addition, the systems with 1400 or more 
particles display a negative minimum. The 700-particle 
system does not show a minimum in D pp (z) indicating 
that the perpendicular energy dissipation in this system 
occurs through particle- wall collisions. Figure [7J\ shows 
that as the number of particles increases, the maximum 
in T> pp (z) increases and the minimum decreases. There 
is also a small shift in the height at which the maximum 
is observed. The shift in the positions of the extrema 
to smaller heights is due to the successive cooling and 
compression of the bed as the number of particles is in- 
creased. 

The dissipation in the parallel direction is always nega- 
tive and exhibits a minimum for all the systems as seen in 
Figure [7)3. The position of the minimum decreases from 
around z = 1.0 to z = 0.5 as the number of particles in- 
creases. The heights of these minima coincide with those 
observed for the maxima in the perpendicular direction. 
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FIG. 6: Contour plots of (A) T>p P (r, z), (B) T>p P (r, z), and (C) packing fraction as a function of height and radial position. The 
contours designate a change of 4.57 in graphs A and B. Simulation conditions are the same as those for Figure [5] 

It should also be noted that the maximum loss in Vpp(z) positive for the two largest velocities, 
is greater than the maximum gain in T> pp (z), as expected 
for inelastic particles. 



D. Influence of the base velocity on local 
dissipation 

We also examined the effect of varying the base ve- 
locity on the dissipation. Systems consisting of 700 
to 4200 particles were simulated with base velocities of 
0.786 < Vb < 15.7. Figure |H] displays the dissipation 
in the perpendicular and parallel directions for a system 
with 2100 particles at each base velocity. An increase in 
the base velocity from 0.786 to 3.14 results in an increase 
in the positive maximum of T> pp (z). There is a corre- 
sponding decrease in the minimum observed in T>p p (z). 
Furthermore, the maximum in the perpendicular direc- 
tion shifts to a larger height as the velocity is increased. 
The shifts in the height are related to the expansion of 
the bed. The velocities of Vb — 7.86 and 15.7 exhibit a 
decrease in the magnitude of the extrema for the perpen- 
dicular and parallel directions. The decrease does not 
translate to a decrease in energy transfer, however, since 
the curves broaden noticeably for the larger velocities. 
The energy transfer is no longer confined to a small re- 
gion of heights, but now extends over a larger region. It 
is also interesting that the values of V p (z) are always 



E. Influence of the particle-particle restitution 
coefficient on local dissipation 

We also collected the values of T>l P (z) and V pp (z) for 
systems in which the particle-particle restitution coef- 
ficient (c) was changed (Figure 13 ■ We varied c from 
0.8 to 0.99 for a system with 2100 particles at a con- 
stant base velocity of Vb = 1.57. The dissipation in the 
perpendicular direction (Figure reaches the highest 
maximum and lowest minimum for the system with the 
restitution coefficient of c = 0.8. These extrema decrease 
in magnitude as the restitution coefficient increases, and 
for c = 0.99 there is no observable minimum in D pp (z). 
Unlike the previous studies of the velocity, these curves 
do not broaden. The peak does not shift and there is only 
a small decrease in height and increase in breadth. The 
dissipation in the parallel direction (Figure[5}3) is always 
negative, as seen previously. The minima become shal- 
lower as the restitution coefficient increases. It should be 
noted that for the nearly elastic system (c = 0.99) the 

magnitude of the minimum in T> pp {z) is nearly equal to 
the magnitude of the maximum in T> pp (z) as would be 
expected for an elastic system. 
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FIG. 7: The dissipation in the (A) perpendicular and (B) par- 
allel directions as a function of height for values of N ranging 
from 700 to 4200 particles (1 to 6 layers of particles). Other 
simulation conditions are R = 14.5, A = 0.348, v — 1.13, 
c = 0.91, and c w = 0.68. 



F. Influence of the particle-wall restitution 
coefficient on local dissipation 

The energy balance, Equation 1211 implies that the en- 
ergy dissipated at the wall will influence Pi, (z). We in- 
vestigated this relationship by simulating systems with 
varying values of the particle-wall restitution coeffi- 
cient. The results for a system consisting of 2100 par- 
ticles and particle-wall restitution coefficients of c w = 
0.5, 0.68, 0.8, 0.99 are presented below. 

The values of Di, (z) and T>\ v (z) as a function of height 
are presented in Figure 1101 The dissipation in the per- 
pendicular direction fFigurelTUlV) does not vary greatly, 
but there is a noticeable trend. The maximum observed 
at small heights becomes less positive and the minimum 
becomes more negative as c w approaches unity. This in- 
dicates that more energy is lost through particle-particle 
collisions as the particle-wall restitution coefficient ap- 
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z 

FIG. 8: The dissipation in the (A) perpendicular and (B) 
parallel directions for systems in which the velocity of the 
base (vb)was varied from 0.786 to 15.7. The simulation 
conditions were N = 2100, R = 14.5, A = 0.348, v = 
0.565, 1.13, 2.26, 5.65, 11.3, c = 0.91, and c w = 0.68. 



proaches unity. The minima in T> pp (z), as seen in Figure 
HOB, indicate that there is less energy loss in the par- 
allel direction as the particle-wall restitution coefficient 
approaches unity. 



IV. CONCLUSIONS 

In this paper we have examined in detail the en- 
ergy transfer between the direction parallel to the en- 
ergy source and the perpendicular directions in a driven 
granular system. Energy balances show that the energy 
changes resulting from particle-particle collisions in these 
directions D pp and D pp are always negative and posi- 
tive, respectively, for inelastic systems. We introduced 
the fractional energy transfer, [3 = —D pp /Dp P , which 
in the steady state equals the fraction of energy input 
at the base that is dissipated at the wall. We further 
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FIG. 9: The dissipation in the (A) perpendicular and (B) 
parallel directions as a function of height for systems in which 
the particle-particle restitution coefficient c was varied from 
0.8 to 0.99. The simulation conditions were N = 2100, R = 
14.5, A = 0.348, v = 1.13,and c w = 0.68. 



examined the quantitative dependence of these quanti- 
ties on the number of particles, the base velocity, the 
particle-particle restitution coefficient, and the particle- 
wall restitution coefficient. We then examined the local 
dissipations, T>p P (r, z) and T>p p (r, z). While the former is 
negative throughout the system, the latter changes sign. 

When an energy source is present in an experimen- 
tal system, it is always anisotropic to some degree, i.e. 
energy is supplied preferentially in some direction(s). 
Therefore the differences in the parallel and perpendic- 



ular directions observed here must always be present in 
real systems. Systems with lower symmetry than the one 
studied here, e.g. a rectangular prism, are expected, in 
addition, to exhibit differences in each non-symmetry- 
related perpendicular direction. 
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FIG. 10: The dissipation in the (A) perpendicular (B) par- 
allel directions as a function of height for systems in which 
the particle-wall restitution coefficient c w was varied. The 
simulation conditions were N = 2100, R = 14.5, A = 0.348, 
v = 1.13, and c = 0.91. 
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